Stability of frozen waves in the Modified Cahn— Hilliard model 
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We examine the existence and stability of frozen waves in diblock copolymers with local conserva- 
tion of the order parameter, which are described by the modified Cahn-Hilliard model. It is shown 
that a range of stable waves exists and each can emerge from a 'general' initial condition (not only 
the one with the lowest density of free energy). We discuss the implications of these results for the 
use of block copolymers in templating nanostructures. 

PACS numbers: 64.70.km, 81.07.-b, 82.35.Jk 



I. INTRODUCTION 

The Cahn-Hilliard equation [l| is often used to model 
microstructures arising from spinodal decomposition in, 
say, polymer mixtures. One of the simplest systems ex- 
hibiting this behavior would be a mixture of two poly- 
mers made from monomers, A and B, with distinct chem- 
ical properties - e.g., if A is hydrophilic whereas B is hy- 
drophobic. In such cases, a monomer unit is attracted to 
units of the same type, while being repelled by the other 
type, implying that the most energetically favorable state 
is the one where A and B units are fully segregated. Such 
a tendency is indeed confirmed by numerical modelling of 
the Cahn-Hilliard equation [2| and is also in agreement 
with stability analysis of similar models Q . 

One feature of the Cahn-Hilliard model is that the or- 
der parameter is conserved globally (reflecting the mass 
conservation law). The standard model, however, can be 
modified for microstructures where the order parameter 
is conserved locally Q. The modified model applies, for 
example, if chains of the A and B monomers are parts 
of the same polymer molecule, known as a 'block copoly- 
mer' [5|j, in which case they can never be separated by a 
distance larger than the size of a single molecule. 

Systems with locally conserved order parameter are 
of particular interest in nanotechnology. In particular, 
block copolymers are used to template nanopatterns at 
surfaces, on scales that are too small for traditional top- 
down photolithography 6]. Such patterns have to be 
'directed' using chemical pre-patterning or topography, 
which is known as graphoepitaxy and can provide ex- 
cellent pattern registry [7j. In its simplest form, conven- 
tional lithographic techniques are used to create trenches 
in a silicon wafer - then the trenches are filled with block 
copolymer which orders into lamellae parallel to the side- 
wall on annealing [8|]. Finally a selected monomer is 
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chemically etched away and the remaining polymer used 
as an etch mask to facilitate pattern transfer to the sub- 
strate, creating nanowires on a scale too fine to be man- 
ufactured by conventional techniques [H, H3] • 

The lamellae used to template the nanowires corre- 
spond to frozen waves (i.e. periodic time-independent 
solutions) of the one-dimensional version of the modi- 
fied Cahn-Hilliard equation. It is particularly important 
whether these solutions are unique or perhaps multiple 
stable solutions exist, as the latter would impede one's 
control over the manufacturing process. 

The present paper answers the above question by ex- 
amining the existence and stability of frozen waves in the 
modified Cahn-Hilliard equation. In Sect. HH we shall 
formulate the problem mathematically. In Sect. IIII[ the 
existence of frozen-wave solutions will be discussed. In 
Sect. IIV1 we shall present the results of a stability anal- 
ysis of frozen waves. 



II. FORMULATION 

Consider a one-dimensional diblock polymer, with the 
characteristic thickness I of the A/B interface and mobil- 
ity M (the latter characterizes the diffusion of the order 
parameter tfi). Using I and l 2 /M to non-dimensionalize 
the spatial coordinate x and time t respectively, we can 
write the one-dimensional version of the modified Cahn- 
Hilliard equation (MCHE) in the form 



+ a<p = 0, 



(1) 



where a determines the ratio of the characteristic size of 
the region over which the order parameter is conserved 
to I. 

As shown in Ref. [ll| , the MCHE admits frozen waves 
only if 

0<a< i, 

whereas the wavelength (spatial period) A must satisfy 
2tt . 2tt 



< A < 



^I(1 + V T=45) yf\ (1 - VT=ttj 



(2) 



(see Fig. 1). Ref. [TT| also computed examples of frozen 
waves and the energy density E as a function of a frozen 
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FIG. 1: The existence region of frozen waves (those described 
by one-wave solutions) on the (a, A)-plane (a is the parame- 
ter in the modified Cahn-Hilliard equation ([T]), A is the wave- 
length). The boundaries of the region are given by 

wave's length A - which turned out to have a minimum at 
a certain A = Ao- Given that the energy cannot grow and 
is, thus, a Lyapunov functional, a conclusion was drawn 
that Ao is stable. 

However, the fact that Ao minimizes the energy of 
frozen waves means only that the corresponding wave 
is stable with respect to perturbations of its length, but 
not necessarily to arbitrary perturbations (for which the 
'general' second variation of E needs to be examined). 
On the other hand, stability rarely occurs for a single 
value of a parameter - hence, if Ao is indeed stable, it is 
likely to be one of a range of stable wavelengths. 

In what follows, we shall clarify the above issues by 
examining the stability of frozen waves through the stan- 
dard linear analysis, not involving energy arguments. 

To do so, we introduce frozen waves as time- 
independent solutions, <j> = 4>{x), for which Eq. ([lj yields 

(4>-4> 3 +$ xx ) xx + a0 = O. (3) 

Together with the condition of spatial periodicity, 

4>{x + A) = $(x), (4) 

Eq. ([3]) determines <fi(x). 

To examine the stability of a frozen wave 4>(x), assume 

< l> = $(x) + $(t,x) ) (5) 

where <fi describes a small disturbance. Substituting ([5]) 
into Eq. ([1]) and linearizing it, we obtain 

4> t + U ~ U 2 4> + kx) +a4> = 0. (6) 

V / xx 

We confine ourselves to disturbances with exponential 
dependence on t (which are usually a reliable indicator 
of stability in general), 

4>(x,t)=^(x)e st , (7) 



where s is the disturbance's growth/decay rate. Substi- 
tution of into © yields 

tnl>+(i/>-34?il> + il>«x) xai + anl> = 0. (8) 

Unlike the base wave </>, the disturbance iji does not have 
to be periodic; it is sufficient that the latter is bounded 
at infinity. Given that ip is determined by an ordinary 
differential equation with periodic coefficients [Eq. ©], 
the assumption of boundedness amounts to the standard 
Floquet condition, 

1>(x + A) = ifi(x) e w , (9) 

where 9 is a real constant. Physically, condition ([9]) im- 
plies that, if the disturbance propagates by one wave- 
length of the base solution, the disturbance's amplitude 
remains the same, whereas its phase may change by a 
value of 9. 

Eqs. ([H])-® form an eigenvalue problem, where s and 
-0 are the eigenvalue and the eigenfunction. If, for some 
values of the phase shift 9, one or more eigenvalues exist 
such that Re s > 0, the corresponding base wave 4>{x) is 
unstable. 



III. FROZEN WAVE SOLUTIONS 

It turns out that a lot of physically meaningful infor- 
mation about frozen waves can be obtained in the limit 
of weak nonlincarity, i.e. under the assumption 

\4>\ < l. 

In order to understand qualitatively what to expect in 
this case, one can simply omit the nonlinear term in Eq. 
@ and seek a solution of the resulting linear equation in 
the form 

4> = e cos {kx + p) , (10) 

where the wave's amplitude e and phase p are arbitrary, 
whereas the wavenumber k satisfies 

-fc 2 (l-fc 2 ) +a = Q. (11) 

Assuming k > and recalling the relationship between 
the wavelength and the wavevector, 



one can see that, in the linear approximation, only two 
wavelengths are allowed - coincidentally, the same values 
which represent the upper and lower boundaries of the 
existence region shown in Fig. 1). Under the weakly 
nonlinear approximation, in turn, one should expect the 
wavelength to be close, but not necessarily equal, to one 
of the above two values, with a deviation from them pro- 
portional to some degree of the wave's amplitude e. So- 
lutions similar to (fit))) are the ones computed in Ref . [Tljj ; 
they will be referred to as 'one-wave solutions'. 
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FIG. 2: A schematic illustrating linear two-wave solutions. 
The solid curve [determined by (|ll[) - (|12[l ] corresponds to lin- 
ear frozen waves. The lengths Aj.,2 of the two 'composite' 
waves are shown by dots and are connected by dashed lines. 
Two-wave solutions exist only for those values of a for which 
A2/A1 is a rational number (presented in the figure above the 
corresponding A2). The wavelength of a two- wave solution as 
a whole equals the lowest common multiple of Ai and A2. 



To understand the physical meaning of two-wave solu- 
tions, seek a solution of the linearized version of Eq. (|3]) 
in the form 



£1 cos (kix + pi) + e 2 cos (k 2 x + p 2 ) 



(13) 



where £1,2 and pi.2 are arbitrary and ki t 2 are the roots 
of Eq. (JTTJ> such that 



< ki < k 2 . 



(14) 



Physically, solution (|13p represents a superposition of two 
periodic waves of zero frequency, with wavenumbers k\ y2 
and phases pi <2 - Note that (fT3)) is periodic only if k 2 jk\ 
(and, hence, A2/A1) are both rational numbers, which 
occurs only for some values of a (several examples of 
such are illustrated in Fig. 2). 

The above linear analysis, however, leaves several im- 
portant questions unanswered. On the (A, a)-plane, for 
example, one-wave solutions seem to exist near curve 
(fTT|) -(fT2 |) . whereas Ref. [ll[ found frozen waves only 
inside this curve, and not outside (see Fig. 1). This 
discrepancy - as well as the question of existence of two- 
wave solutions - will be clarified in Sects. IIII AI IhTTJI 
The case of strong nonlinearity will be examined numer- 
ically for both types of frozen waves in Sect. IIII CI 



A. Asymptotic results: one-wave solutions 

It is convenient to change the spatial coordinate x to 
£ = kx, 

[where k is, again, the wavenumber determined by (|12[) ]. 
In terms of £, Eqs. ©-((U) become 



+ fc 2 0£f)^ + a4> = 0, 



(15) 



4>{x + 2tt) = cj){x). (16) 
We shall seek a solution as a series of the form 

= £(V O) +£ 2 ^ 2) + •■■), (17) 

and also expand the wavenumber k, 

k 2 = + e 2^(2) + . . . _ ( 18) 

To leading order, Eqs. (|15l) (fTTJ)) reduce to 

KM + + atf® = 0, (19) 

4>^{x + 2it) = ftV(x). 

We seek a solution in the form 

(o) = Acos«+p), 

where A and p are real constants and n > is an integer. 
It is sufficient to examine the case n = 1, as n > 2 
corresponds to re-defining the solution's spatial period 
by including more than one wavelengths in it, without 
changing anything physically. We also let A = 1 (as the 
wave's physical amplitude still remains arbitrary due to 
the arbitrariness of e) and p = (as a phase constant 
can always be included later). Thus, the leading-order 
solution becomes 



(o) = cos£- 
Substitution of (J2DJ) into ([TSJ yields 

= i (1 ± VT= 15) . 

In the next-to- leading order, Eq. (|15[) yields 



(20) 



(21) 



+ K (0) (IV) _ 0(0)3 + K W$W + 



€€ 

a^ (2) = 0, 



which, upon substitution of (j2"0")) , becomes 



S.S. 



4 

9K(° 



cos3x, (22) 



This equation can have a 27r-periodic solution only if 
the term involving cos£ on the right-hand side vanishes, 
which implies 



4 (l - 2K (°)) ' 



(23) 



4 



after which ((22|) yields 

4> {) = 77^ 77m cos 3f . (24) 

V 36K (°) (1 - 9#( )) - 4a ^ V ; 

Recalling that K is related to the wavenumber k by (fTS)) 
- and fc is, in turn, related to the wavelength A by (fl"2|) . 
one can use expressions (|2"Tj) and f[23]) to obtain 

A 1|2 = - 2 ^ (l± ^ ) + Q(e% (25) 

where i and 2 correspond to + and — respectively. For 
e = (i.e., in the linear limit), Ai and A2 are the lower 
and upper boundaries of the existence interval in Fig. 1, 
whereas an increase in e causes an increase in Ai and a 
decrease in A2. This means that one- wave solutions do 
not exist above (A 2 ) £=0 and below (Ai) - which agrees 
with the conclusions of Ref . [ll[ . 

Observe that expansion (|25l) fails in the limit a —> ^. 
This case should be examined separately, by assuming 

a = \+e 2 a^. (26 ) 

An expansion similar to the one above yields (technical 
details omitted) 

k 2 = § ± ^-e 2 ^ 2 ) - f e 2 + 0(e 2 ), (27) 

from which A can be readily found. The corresponding 
expansion of the solution is (technical details omitted) 

= £[cosC- T fg£ 2 cos3C + O(£ 4 )] . (28) 

Formulae (|27[) -- (j28|) describe the one- wave solution near 
the 'tip' (a — > j, A — > 2^/2tt) of the existence region 
shown in Fig. 1. Formula (f25j). in turn, describes the 
solution near the boundary of the region, but not too 
close to its tip. 

Finally, note that substitution of expression (|2"Tj) for 
K (o) 

into expression ((24)) shows that <f>^ becomes infinite 

for 

a = TUo' ^ = 2VT07T. 

This is another example where the straightforward ex- 
pansion for the one- wave solution is inapplicable (in ad- 
dition to the 'tip point', a — j, X — 2y/2n). The un- 
bounded growth of the first-order term involving cos 3£ 



suggest that, in this case, cf>^ should include both cos£ 
and cos 3£ - and not just the former term, as solution 
pO)) does. 

This case will be examined in the next section. 



B. Asymptotic results: two-wave solutions 

As mentioned above, two-wave solutions exist near 
those values of a for which the ratio of the wavenum- 
bers of the two individual waves is a rational number. 
Let a(°' be one of such values, with being a small 
deviation from it, 

a = a^+s 2 a^. (29) 

Substitution of this expression, together with expansions 
fT7 ]l -([T8" |) . in Eq. ([15]) yields, to leading order, 

K<® (0 (o) + if (0) 4?) a + « (0) ^ (0) = °- ( 3 °) 

We assume that the two waves of which a two-wave so- 
lution consists have wavenumbers K^°'rii and K^n2, 
where 711,2 are integers (without loss of generality, they 
can be assumed to be coprime and such that n\ > n-i > 
0). Accordingly, we shall seek a solution of Eq. (|3U)) in 
the form 

= a cos(nix + 6) + cos n%x (31) 

(the fact that the second wave's amplitude is unity and 
its phase is zero does not reduce generality). Substituting 
flEE) into ([5D]). we obtain 

a(°>= nlnl 2 , A-(°) = ^i-_ (32) 

Recalling that K is related to the wavelength A, one 
can see that these equalities determine the points in the 
(a, A)-plane near which two- wave weakly-nonlinear solu- 
tions are localized (they are illustrated in Fig. 2). 
To the next-to-leading order, ([15)) yields 
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K<® (^ (2) + ^ (O) 0g ) ) 5? + « (O) (2) = cos(n lX + 0) + ^ [3cos(nia! + 0) + cos3(thx + 0)] 

3a 2 3(2 
+ — {cos [(2ni - n 2 ) x + 20] + cos [(2ni + n 2 ) x + 20]} + — {cos [(m - 2ti 2 ) a; + 0] + cos + 2n 2 ) x + 9]} 



H cos n 2 x H — (3 cos n 2 x + cos 3n 2 x) 

2 4 



£T (2) n 2 AjA" (0) n 2 - l) + a (2) cos(nix + 0) 



K^n 2 [2K^nl-l)+a^ 



cos n 2 x. 



This equation has a periodic solution for </>( 2 ) only if the the specifics depend on whether (m, n 2 ) = (3, 1) or not. 
terms involving cos(nix + 0) and cosn 2 a; cancel out, but In the former case, straightforward calculations yield 



(2) _ 3 (7a 3 - a 2 + 17a + 3) (2) _ 9 (9a 3 + 3a 2 + 9a + l) 



340a 



420a 



if (ni,n 2 ) = (3, 1) , 



whereas the latter case yields 

3 [{n\ - 2n|) a 2 + 2n\ - n 2 ] 



4 (n 2 + n\) {n\ - n\) 



(2) g»fea (« 2 + 1) - f 
a> ' = — -^r- it n\ T 3n 2 . 



4 (n\ + n\) 2 



r 



Now, using expressions (fT7 ^ — ([TH ]) . {29]) , (fT2|) . (|3T ]i -([32" l) and 
to relate 4>^ 2 \ -ft^ 2 -*, and a^ 2 ' to the 'physical' quantities, 
we can summarize the two-wave solution in the form 



■ e [a cos(nia; + 0) + cos n 2 x] + 0(e 3 ), (33) 



9 

100 



,9 (9a 3 + 3a 2 + 9a+ 1) 



A = 2y/Wn + e- 



+ ' e4 »- 

,37rv / T0(7a 3 - a 2 + 17a + 3) 
34^ 



+ 0(e 4 ), 



if (n 1 ,n 2 ) = (3, 1) , 



(34) 



(n 2 + n 2 ) 2 



_ £2 9n 2 n 2 (a 2 + l )+o(£4 ^ 



4 (n 2 + n|) 



A = 2-K^/n\ + n| + e 



2 3tt [(n 2 - 2n 2 ) a 2 + 2n 2 - n 2 ] ^Jn\ + n 2 + Q 



> if m ^ 3n 2 . 



4 (nf - n 2 ) 



(35) 



Expressions (f34|) and (|35|) can be viewed as parametric 
representations (with e and a being the parameters) of 
the existence region of two- wave solutions with ni j2 , on 
the (a, A)-plane. It can be shown that, for all ni )2 , this 
region is a 'semi-infinite sector' (see an example in Fig. 



3). However, since we assumed weak nonlinearity, this 
conclusion can only be trusted near the vertices of the 
sectors. Effectively, we have found the tangent lines to 
the boundaries of the 'true' region of the existence inter- 
val. 
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FIG. 3: Regions of existence (bounded by the dotted line) of 
two-wave solutions with A2/A1 = 2:1. Thin solid lines show 
the boundaries of the existence region calculated under the 
assumption of weak nonlinearity. The thick solid line shows 
the existence region for one- wave solutions (as in Fig. 1). 
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100 



6 (9a 3 + 3a 2 + 9a + l) 



A - 2VT0tt 27%/TOtt (7a 3 



+ 17a + 3) 



+ 0(e 2 ), if (ni,n 2 ) = (3,l), 



(36) 



or 



(n? 



§n\n\ (a 2 + l) (n\ - n|) 



A - + n'i 3ir {n\ + n 2 2 f /2 [{n\ - 2n 2 ) a 2 + 2n\ 



+ 0(e 2 ) if m ^ 3n 2 . 



(37) 



For given m,a and (a, A), (|36|) and (|37|) can be treated 
as equations for a [it is, essentially, the ratio of the am- 
plitudes of the waves which constitute the two-wave so- 
lution - see (|33p]. Observe that, (|3"7| admits two roots 
for a (with equal magnitudes and opposite signs), which 
means that two-wave solutions with n\ 7^ 3ri2 exist in 
pairs. In some cases these solutions can be obtained from 
each other by shifting x — > x + const (for the 2:1 case, 
for example, const = tt) — but in other cases, the solu- 
tions with positive and negative a seem to be genuinely 
different. 



Numerical results 



In this section, we shall present examples of strongly- 
nonlinear one- and two- wave solutions and the existence 
region for the latter. Two numerical methods have been 
used: the shooting method (which turned out to be insuf- 
ficient for large wavelengths) and the method of Newton 
relaxation (which worked marginally better). 



For the case (711,712) = (3,1), in turn, it follows from 
(|3"T|) that three roots exist for a. One of the three, how- 
ever, corresponds to the one-wave solution. Indeed, recall 
that the expansion derived for those failed near a = , 
A = 2 V / T07r, i.e. precisely where the asymptotic theory 
for too-wave solutions predicts existence of those with 
(711,712) = (3,1). In fact, one- and two-wave solutions 
cannot be distinguished in this region, as, in both cases, 
the coefficients of cos3x and cos a; are of the same order. 



Figs. 4a and 4b show examples of increasingly nonlin- 
ear one- and two-wave solutions respectively, for a fixed 
wavelength A and decreasing a. In Fig. 4a, observe the 
increase of the wave's amplitude as a is moving away 
from the boundary of the existence region (for this value 
of A, the boundary is located at a m 0.131). Fig. 4b, 
in turn, illustrates the fact that the margins of the ex- 
istence region for two-wave solutions correspond to the 
cases where the amplitude of one of the two waves van- 
ishes (which is how they bifurcate from one-wave solu- 
tions). 
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FIG. 4: Examples of frozen waves with spatial period A = 
16: (a) one- wave solutions with a = 0.005, 0.050, 0.100 
(curves 1, 2, 3 respectively); (b) two-wave solutions with 
a = 0.07, 0.09, 0.11 (curves 1, 2, 3 respectively). 

IV. THE STABILITY OF FROZEN WAVES 
A. Asymptotic results 

We shall first examine the stability of frozen waves 
asymptotically, under the same assumption of weak non- 
linearity used to find the frozen waves themselves. 

We shall start by re-writing the stability problem @- 
([9]) in terms of £ = kx, which yields 

sip + k 2 (tp - 34> 2 ip + fcVft)^ + atp = 0, (38) 

<K£ + 2tt) =V(£) e ie . (39) 

It turns out that all weakly nonlinear solutions, one- and 
two-wave alike, are unstable for all a except a — > 7 (i.e. 
except the 'tip' of the existence interval for one-wave so- 
lutions shown in Fig. 1). We emphasize that this conclu- 
sion does not apply to strongly nonlinear waves (which 
can be either stable or unstable for any a - see the next 
section). 

Thus, assuming that a = j + 0(e 2 ), we represent a 
by expression (I26|) ; the corresponding base solution, 0, is 
represented by ([2T | -(|28" ]l . 



A straightforward analysis shows that, if the phase 
shift 9 is order-one, all solutions of the eigenvalue prob- 
lem ([3"8]) - (j39|) . with the parameters determined by (|26])- 
(f2"5|) . are stable. Thus, instability can occur only for small 
9, which can be conveniently accounted for by letting 

9 = 2neq, (40) 
where q is an order-one constant. Note that the smallness 
of the phase shift 9 implies that the instability occurs at 
wavelengths close to that of the base solution. 
It is also convenient to introduce 

i> new = i>e-' leqx ■ (41) 

Using (|gQ" ]> -([3I |) to replace 9 and V in Eqs. (|38 ]l -(|3"9" l) 
with q and i/Wu>, we obtain (subscript new omitted): 

sip + k 2 (ip + k 2 ip(.(. - 3e 2 4) 2 i/j)^ 
+ ieqk 2 (2V»{ + ik 2 ip m ) - e 2 q 2 k 2 + 6k 2 ip)^ + 0(e 3 ) 

+ aip = 0, (42) 

<K£ + 27r) = V(0, (43) 

where the specific form of the terms 0(e 3 ) will not be 
needed. 

It can be demonstrated that no instability occurs if 
s = 0(1) - hence, only small s need to be examined. 
We assume (and shall eventually justify by obtaining a 
consistent asymptotic expansion) that 

s = e V 2 ) + ... , 

while the eigenfunction is 

ij) = i' {0) + eV (2) + • • • • 
To leading order, Eqs. ([32" ]) -(|35 1) . l(25 )l -([25 )l reduce to 

i(^ + ^ ) ) ee + i^ = o, 

^ 0) (C + 27r)=V (0) (0, 

which yield 

?A (0) = Ssiax + Ccosx, (44) 

where S and C are undetermined constants. 

In the next-to- leading order, Eqs. flU])-®, (j2"o) - ([2gj) 
yield (after straightforward algebra) 
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(2) 1 . (2) 



[5 sin x + C cos x + ^5* (3 sin 3x — sin x) + \C (3 cos 3x + cos x)] 



s (2) + g 2 _ 3 j (5 s i nx + Ccosx) - 2W-a( 2 ) - | (5 cos x - Csinx) = 0, 



lf»< 2 >(£ + 27r) = ^ w (0. 



r 



This boundary- value problem has a solution for -0^ 2 ' only 
if 



S + («W + c? 2 - §) S + 2iqy/-aW - §C = 0, 



s (2) + ^_3 )c ,_ 2i / ( 2 )_3 S=0) 



which, in turn, has a solution for 5 and C only if 



,(2)2. 



(2g 2 + |) S (2) 



g 2 + 3_ 4 _ a (2)_3 



= 0. 



This equation determines the eigenvalue s^ 2 K It can be 
readily shown that s 1 - 2 ' is stable (i.e. Res' 2 - 1 < 0) for all 
q only if 



9_ 

16- 



Finally, using (|26[) - (|27|) to express in terms of the 
'physical' parameters a and A, we obtain the follow- 
ing stability criterion for a frozen wave with parameters 
(a, A): 



\-2V2tt) <§7r 2 (i 



(45) 



Observe that, for any a < \, there exists an interval of 
stable wavelengths, not just a single value of A. 

As mentioned before, condition (|45|) applies only if a 
is close to \. In the next subsection, it will be extended 
numerically to arbitrary values of a. 



B. Numerical results 

The eigenvalue problem ([8])- ((9]) was solved numerically 
for if)(x) and s with the base wave 4>{x) computed using 
problem 

The general features of the dispersion relation [the 
dependence s(0)] of the eigenvalue problem (j8|)-(j9]) is 
described in the Appendix, whereas here we shall only 
present the stability diagram on the (a, A) plane - see in 
Fig. 5. Evidently, for all values of a, an interval of A ex- 
ists where one- wave solution are stable - which confirms 
and extends the asymptotic (a j) stability criterion 
(|4"S"|) . We have not found any stable two- wave solutions, 




FIG. 5: The stability region of one-wave solutions on the 
(a, A)-plane (shown by the dashed line). The thick solid line 
shows the existence region of one- wave solutions (as in Fig. 1). 
The upper panel represents a blow-up of the shaded region of 
the lower panel [the thin solid line shows the stability region's 
asymptotic boundary, (|45|)]. 



which suggests they are either unstable or perhaps their 
regions of stability are small and difficult to locate. 

To illustrate that any solution from the range of stable 
frozen waves (not necessarily the wave with the minimum 
free-energy density) can emerge from a 'general' initial 
condition, we have carried out the following numerical 
experiment. The time-dependent MCHE ([TJ was simu- 
lated using finite differences with a fully implicit back- 
wards Euler method, and the results of the simulations 
presented below are for a = 0.1. 

In this case, the energy minimizing wavelength is 
Aq = 11.31 (calculated from the interpolation formula 
of Ref. [HI), with the corresponding frozen- wave solu- 
tion denoted by cf>o (x). According to our analysis, how- 
ever, a solution <f>\ (x) with the commensurate wavelength 
Ai = 14.13 (A /Ai = 4/5) should also be stable. 

To verify this, Eq. ([!} was simulated in a domain of 
size L = 5Ao = 4Ai which accommodates both solutions. 
The initial condition was chosen as a 'mixture' of the 
frozen waves <po (x) and <j)\ (x), i.e. 



(46) 



where (3 € (0, 1) is the 'mixing ratio'. The timestep was 
0.1 and 400 gridpoints per period were used, and it has 
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FIG. 6: Evolution of perturbed steady states within the sta- 
ble region at a = 0.1. The traces show the limiting (t — y oo) 
solution initialized by (|46p for various j3. The traces are in- 
cremented in steps of 0.5, beginning with zero increment for 
/3 = (the lowest trace) and ending for j3 — 1 (the highest 
trace). The curves are marked with the corresponding values 
of 13. 

been verified that the results were mesh and timestep 
independent. 

If 0o (x) was the only stable solution, the system would 
evolve towards <pQ for all Our simulations nevertheless 
show that, for (3 = 0.2, the system evolves back to <f>i 
(see Fig. 6), which confirms our conclusion about the 
existence of multiple stable states. 
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e 

FIG. 7: The dispersion relation s(9) for the eigenvalue prob- 
lem ©-(HJ) with a = 0.24, A = 8.886 (solid line). The dotted 
line shows the limiting dispersion curve (IA2() . Note the differ- 
ence in the vertical axes' scales of the upper and lower panels. 

selected rather than a true time independent solution to 
the modified Cahn-Hilliard equation. 

Finally, it would be interesting to extend the present 
results to steady states with two spatial dimensions, sim- 
ilar to those found in Ref. [l2| for an equation similar to 
the two-dimensional MCHE (but with a slightly different 
nonlinearity) . 



V. CONCLUDING REMARKS 

The main result of the present paper is illustrated in 
Fig. 5, which shows the stability region of (one-wave) 
frozen solutions of the modified Cahn-Hilliard equation 
(HJ) . We have also found a new class of frozen waves - the 
'two- wave solutions', but these seem to be unstable and, 
thus, less important than the usual, one-wave type. 

We have also made a more general - and potentially 
more important - conclusion regarding the energy ap- 
proach to studies of stability. If a family of solutions 
exists and one of them minimizes the energy functional, 
this does not necessarily mean that all other solutions are 
unstable. Furthermore, the stability of the minimizer so- 
lution cannot be guaranteed either: even though it is 
stable with respect to the perturbation of 'shifting along 
the family of solutions', another perturbation can still 
destabilize it. 

Physically, our results imply that when lamellar mi- 
crostructures of block copolymers are used to template 
nanowires, one must ensure that only the desired state 
is created. This may become more of a critical concern 
as larger numbers of nanowires are to be created within 
a single trench. In practice, some control over this can 
be exerted via the annealing schedule. It should also be 
noted that a kinetically stable quenched state may be 
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Appendix A: The structure of the dispersion 
relation and the eigenfunctions of the eigenvalue 
problem fl!])-® 

First of all, it can be readily shown (and has been 
confirmed numerically) that the eigenvalue s of problem 
©-© is real. 

Observe also that, for the boundary points of the ex- 
istence region (see Fig. 1), frozen waves have zero am- 
plitude, i.e. = 0. In this case, the eigenfunction of 
problem ([S])-© can be readily found, 

iOx 

XP ~' ' 

as well as the dispersion relation (i.e. the dependence of 
the eigenvalue s on the phase shift 8), 
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2tt (provided they are extended to negative 8). 

As a result, the only 'original' part of the upper branch 
is the segment 8 £ [0,2tt], whereas the 'original' part of 
the lower branch is that for 8 = [2tt, oo). Not surprisingly, 
the unperturbed dispersion relation (IA2|) 'switches' from 
the upper branch to the lower one near the point 8 = 2tt 
(see Fig. 7). 

Fig. 8 shows typical behavior of the eigenfunctions: 
for 9 = 0, tp{x) does not oscillate at all; for 8 = 2tt, it 
oscillates once; for 8 = 4ir, it oscillates twice, etc. 

Fig. 9, in turn, shows the onset of instability brought 
by a change of the period of the base wave. One can 
see that the waves with 8 « 2ir are first to lose stability 
(which agrees with our asymptotic analysis of the case 
aH>±). 

Finally, we mention that two exact solutions were 
found for the eigenvalue problem ((8|)-(j9|): 







for 
for 



0, 
2tt. 



(A3) 
(A4) 



FIG. 8: Examples of eigenfunctions of problem (H}-© with 
a = 0.24, A = 8.886. The eigenfunctions marked with "U" 
and "L" correspond to the upper and lower branches, respec- 
tively, of the dispersion relation illustrated in Fig. 6. 



Observe that, as follows from (|Al[) . the eigenfunction be- 
comes more and more oscillatory with increasing 9. 

In the general case, (i.e. for interior points of the ex- 
istence region of frozen waves, where 4> 7^ 0), one would 
expect that (|A2[) is somehow perturbed, but still keeps its 
structure as a single curve on the (8, s) plane. Further- 
more, the large-# part of the 'general' dispersion relation 
should not differ much from that of (|A2j) - as, for rapidly 
oscillating tp, the term involving in equation ([5]) is neg- 
ligible. 

Our computations show, however, that, no matter how 
small (f> is, the dispersion relation splits into two branches 
- see Fig. 7. 

The following properties of the two branches have been 
observed: 

• The upper branch is periodic with a period of 47T. 

• Both branches are symmetric with respect to 8 = 



In the latter case the disturbance can be found analyt- 
ically, tp — <j), and it corresponds to infinitesimal shift 
of the base wave. The former solution does not seem 
to have an obvious physical meaning (nor does it admit 
an obvious analytical expression for the eigenfunction, as 
equality (|A4j) has been established numerically). 




FIG. 9: The dispersion curves (upper branches) for a = 0.24 
and (1) A = 8.886 (stability for all 8), (2) A = 9.786 (instabil- 
ity for sufficiently small 9). 
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